Prediction of Stable Walking for a Toy That Cannot Stand 
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Previous experiments [M. J. Coleman and A. Ruina, 
Phys. Rev. Lett. 80, 3658 (1998)] showed that a gravity- 
powered toy with no control and which has no statically sta- 
ble near-standing configurations can walk stably. We show 
here that a simple rigid-body statically-unstable mathemati- 
cal model based loosely on the physical toy can predict stable 
limit-cycle walking motions. These calculations add to the 
repertoire of rigid-body mechanism behaviors as well as fur- 
ther implicating passive-dynamics as a possible contributor to 
stability of animal motions. 



Introduction. For walking and other activities, people 
and animals move in complex, yet stable ways. One view is 
that such coordination is the action of neuromuscular con- 
trol constrained by, among other things, the laws of classical 
mechanics. However, one might ask how much of animal co- 
ordination might be understood in purely mechanical terms. 
Likewise, how much versatility of motion is possible with sim- 
ple mechanical devices? This paper concerns one example 
that sheds a little light on these two general questions. 

McGeer's (e.g. §) success with straight-legged, two-dimen- 
sional uncontrolled and gravity powered walking mechanisms 
highlights the possibility of pure mechanics generating coordi- 
nation. McGeer found steady walking solutions (periodic gait 
or limit-cycle motions) that were exponentially stable (asymp- 
totically returned to the periodic motion after small distur- 
bances from that motion) . In his two-dimensional theory, only 
fore-aft stability, and not lateral balance, is in question. In 
his physical implementations side-to-side balance is enforced 
by duplicate side- by-side legs (4 legs total). These machines 
cannot stand fully upright, but can stand with splayed legs, 
possibly contributing to their dynamic stability. 

Extending McGeer's ideas, Coleman and Ruina |^ de- 
scribed an easily reproducible |^ two legged gadget built from 
Tinkertoys® that cannot stand at all, even with both feet on 
the ground, splayed or not, yet seems (slightly) dynamically 
stable while walking. But where the stability of McGeer's ma- 
chines was first predicted with rigid-body modeling, the sta- 
bility of the Tinkertoy® device was not. As noted in [^, this 
system is essentially dissipative (from collisions and possibly 
from ground friction and internal dissipation) and nonholo- 
nomic (the dimension of the accessible configuration space is 
larger than the dimension of the velocity space). Nonholo- 
nomic systems can have asymptotic stability, even when con- 
servative, and nonholonomicity from intermittent foot contact 
might also contribute to stability |^,^]. 

Which properties are needed for asymptotic dynamic sta- 



bility of such a statically unstable system was left unanswered 
by H. Possible key factors include friction of the hinge, play 
in the hinge joint, elastic or inelastic deformation of the struc- 
ture, compliance at the foot contact, and sliding and twisting 
friction at the foot contact. Could a rigid body model without 
these effects explain the stable motion? 

Previous research. The simulation model in con- 
sists of two rigid bodies connected by a frictionless hinge 
(Fig. |l|). The feet are toroidal with principal radii ri and r2. 
The ground allows no relative motion of contacting points, 
no torques at the foot contact, no bounce (restitution e = 0), 
and no tension force at the foot contact (non-sticky floor). 
The lengths, center of mass location, the moment of inertia 
components, and the ground slope are adjustable. After non- 
dimensionalizing with mass m, length / and time \JTfg there 
are 13 free parameters. The working Tinkertoy® was based 
on mildly unstable simulations Q of a simplified model with 
point contact (ri = r2 = 0) and narrow hips (ui = 0). 

Earlier, McGeer ^ studied the same model, allowing r\ > 
and w > ^ but assuming that the principal moments of in- 
ertia aligned with the hip hinge and leg. He found only unsta- 
ble solutions where, also, the swing foot passed below ground. 
Kuo studied a similar model, but disallowed steer ((^) and 
also only found unstable passive gaits. Dankowicz found 
stable solutions for a related kneed computational model. 
That model has wide feet so, like the 2D models, it can stand 
stably with splayed legs. The semi-3D computational model 
of Wisse et al. Q can also stand with splayed legs. 

Methods. Our study was of the system in (Fig. |l|) [^] 
described above, but with hip spacing and disc feet \r\ 7^ 
0,iu 7^ 0, r2 = 0). 

The overall approach is to characterize the solution of the 
rigid-body dynamics equations as a function (map) that takes 
the state of the system just after one step as input and gives 
the state just after the next step as output. A fixed point 
of this map defines a limit cycle. Stability is evaluated by 
the eigenvalues of the matrix describing the linearization of 
this map. If all eigenvalues are less than one in modulus 
the periodic motion is exponentially stable. The map, its 
fixed points and the linearization are all found from numerical 
solutions of the equations of classical rigid-body dynamics. 

The numerical searches were aimed at generating stable 
motion and not at accurately modeling either the existent 
physical toy or humans. We used the toy's approximate pa- 
rameters to seed the optimizations. Special purpose optimal 
control software (see Appendix) was used to reduce the most 
unstable eigenvalue while maintaining periodicity of the gait, 
positive foot clearance and static instability. The resulting so- 
lution was checked and improved with an independent method 
and checked again with another independent simulation. 
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Results. The model of Fig. |l| has asymptotically sta- 
ble limit cycle motions (Fig. |2|), with the foot of the swing 
leg clearing the ground, with Ixx ~ 0.1982, Iyy = 0.0186, 
Izz = 0.1802, IxY = 0.0071, Ixz = -0.0023, Iyz = 0.0573, 
a = 0.0702, Xcm = 0, Fcm = 0.6969, and Zcm = 0.3137, 
W = .3624, and Ri = 0.1236 and R2 = 0. Capital letters 
indicate non-dimensional variables. Tensor components Imn 
and mass positions are in local left-leg coordinates with origin 
at the vertically-standing contact. Note the static instability 
{Zcm > Ri)- The largest eigenvalue modulus of the single- 
step map Jacobian is 0.8391560, safely below 1. 

Discussion We are claiming a qualitative theoretical 
mechanics result. That is, a system described with the clas- 
sic equations of rigid body mechanics has an exponentially 
stable limit-cycle solution in the neighborhood of a statically 
unstable configuration, but with no fast spinning parts. Al- 
though we have not attempted a mathematical proof, we have 
attempted to do our numerics with enough checks and tests 
to state the result with confidence (see appendix below). 

This solution can be interpreted as bipedal walking, al- 
though not especially anthropomorphic. The base solution is 
exactly repetitive, step after step. That the largest eigenvalue 
is less than 1 in magnitude means that, if the mechanism were 
slightly disrupted from this periodic motion it would asymp- 
totically approach this motion again, over a number of steps. 

Simple numerical probes show, as exponential stability de- 
mands, a small but non-infinitesimal basin of attraction. We 
have not investigated the shape or size of this domain in de- 
tail; we do not know exactly what set of motions eventually 
settle into the periodic motion and thus cannot precisely de- 
scribe what disturbances knock the walker down. However, 
the success of the physical model |^ suggests that the stability 
is robust enough to be physically relevant. 

We do not claim to have an accurate model of the toy in 
[^. Rather we have a simple model that explains the toy's 
qualitative behavior. Accurate quantitative prediction of the 
toy's motions may well depend on physical effects that are not 
in our simple model (various frictions and deformations) . Yet 
unknown is whether the parameters presented here could be 
used as a basis for a better working physical device. More 
generally, we also do not know if more human-like stable 
passive-dynamic designs can be made which are also stati- 
cally unstable. 

Conclusions. The dynamic stability of a statically un- 
stable walking mechanism can be predicted with a model con- 
sisting of two rigid bodies connected by an ideal hinge and 
making intermittent ideal no-slip, no-bounce point contact 
with the ground. We have shown that there is no need to ap- 
peal to hinge-friction, hinge-play, distributed or contact defor- 
mation (elastic or inelastic) , or contact frictional slip in order 
to qualitatively predict the interesting behavior demonstrated 
by the physical model in 

The results further highlight the versatility of simple pas- 
sive strategies for stabilization of coordination. The calcu- 
lation also slightly expands the range of known rigid-body 
phenomenology. 

For videos and reprints about the Tinkertoy® and related 
machines, see www.tam.cornell.edu/~ruina/hplab/pdw.html. 
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Appendix: Numerical Analysis We carried out the 
numerical analysis three different ways. 

The first stable solution (with |crmai| = 0.897) with foot 
clearance was found using the approach developed by Mom- 
baur et al. ||l^ on the basis of the optimal control code Mus- 
cod by Bock et al. and Leineweber jl2j. In the language 
of the discipline, Muscod has been written for general multi- 
phase optimal control problems and is based on a multiple 
shooting state discretization. Multiple shooting splits up the 
original boundary value problem into a number of initial value 
problems enforcing continuity conditions at the transitions 
from one interval to the next. At all multiple shooting points 
Muscod allows the user to impose a number of equality and 
inequality constraints on the parameters and the dynamic 
variables being varied in the optimization. For the tinkertoy 
model described here, this permitted us to ensure periodicity, 
foot clearance during the step, and to keep all state variables 
and parameters within reasonable ranges. Sensitivities of the 
integration end values on each interval, both to variations 
in initial values and to variations in model parameters, are 
efficiently computed by means of internal numerical differen- 
tiation (IND). The basic principle of IND is to use identi- 
cal, but adaptive and error-controlled discretization schemes 
for integration and derivative generation. For use with actu- 
ated and passive gait problems, with implicit state-dependent 
phase switching points and discontinuities in the state vari- 
ables, the original Muscod has been combined with an object 
oriented modeling library that deals with these situations in 
a uniform and complete way. Also added to Muscod were sta- 
bility analysis modules which compute the linearized Poincare 
map of a periodic solution and their eigenvalues, assembling 
information from all multiple shooting intervals and taking 
care of the above mentioned implicit switching point depen- 
dencies. Stable solutions for the tinkertoy model were found 
by varying model parameters and bounds based on coarse grid 
sensitivity information gathered during the previous optimal 
control problem solutions. 

Second, we reproduced and improved the solution above by 
4th order Runge-Kutta integration of the ODE's in Matlab® , 
finding the collision time accurately using Henon's method 
(changing the independent variable near the collision time). 
The impact jump transition is a matrix multiplication. The 
fixed points of the resulting return map were found by numer- 
ical root finding. The fixed point map Jacobian was found 
by a central difference perturbation of the initial state. The 
eigenvalue was reduced from 0.897 to 0.839 using a simulated 
annealing optimization of the maximum eigenvalue modulus. 
For this Matlab solution we did extensive convergence tests 
on both the integration step size and the central difference 
step size. These tests indicate a combined roundoff and trun- 
cation error of about ±10~^ in the largest eigenvalue modulus 
for the parameters given. This maximum eigenvalue estimate 
differs from that generated by Muscod with these parameters 
by 2 X 10"^ 

Finally, the equations of motion were derived indepen- 
dently and simulated again independently in Matlab giving 
agreement to the Matlab solution above to 10~^ for the largest 
eigenvalue modulus. 

For reference, the state of the system just after colli- 
sion of the left foot is, for the parameters given, q* = 
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[0,7/>,6l,t,e,„,<j!>,V,6'.t,e.»] = [0.09868, -0.00925, -0.16016, 
3.43583, -0.13221, -0.01991, 0.47124, -0.39256] with step 
period r* = 1.00711 where ( ) = d( )/dr with r the dimen- 
sionless time. 
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T, non-dimensional time 

FIG. 2. A gait cycle (two steps). In the first half, the 
stance leg is the left leg. (a) The swing leg angle is here mea- 
sured from the slope normal {Os^{t) + 6tt{T) — tt); (b) pos- 
itive swing-foot clearance between collisions; (c) the motion 
has more steer (yaw) than lean (bank). 
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